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The phenomenon of spatial clustering induced by death and reproduction in a population of 
anomalously diffusing individuals is studied analytically. The possibility of social behaviors affecting 
the migration strategies has been taken into exam, in the case anomalous diffusion is produced by 
means of a continuous time random walk (CTRW). In the case of independently diffusing individuals, 
the dynamics appears to coincide with that of (dying and reproducing) Brownian walkers. In the 
strongly social case, the dynamics coincides with that of non-migrating individuals. In both limits, 
the growth rate of the fluctuations becomes independent of the Hurst exponent of the CTRW. The 
social behaviors that arise when transport in a population is induced by a spatial distribution of 
random traps, have been analyzed. 



PACS numbers: 02.50.Ey,05.40.Fb,87.23.Cc, 



I. INTRODUCTION 

A population of Brownian walkers in the presence of 
processes of birth and death, is known to undergo phe- 
nomena of spatial clustering [1, 2]. It was argued that 
this effect could contribute in important way to the spa- 
tial inhomogeneity of the plankton distribution in the 
ocean [3], and could have application in the modeling of 
algal blooms. 

The mechanism is of quite general nature, and can be 
seen as a simple "bosonic" realization of directed per- 
colation [4]. There are several situations, however, in 
which normal diffusion is not an appropriate description 
of individual migration. In the case of plankton, for ex- 
ample, the presence of long lived vortices may lead to 
trapping phenomena and to the possibility of subdiffusivc 
behaviors [5] . Similar coupling of birth-death processes 
and subdiffusion are expected to occur in the spreading 
of radioactive substances in geologic structures [6] . On 
the other hand, the migration strategies of living organ- 
isms, over a large range of scales, is often associated with 
long jumps [7, 8], and this may lead to superdiffusion [9]. 
Again, the generalization from Brownian to anomalously 
diffusing agents has a close counterpart in the anomalous 
directed percolation studied in [10], in which jumps are 
modeled by Levy flights. 

It should be stressed that modelling transport in a ran- 
dom environment as a diffusive process, even in the ab- 
sence of demography, is a non trivial problem (see e.g. 
[11, 12]). The correlations induced by the spatial struc- 
ture of the environment are clearly lost if the individuals 
are treated as independent random walkers. In the case 
of living organisms, these correlations could be seen as a 
form of social behavior, and this leads naturally to the 
question of what is the impact of social behaviors, in the 
broad sense, on transport. 

We have considered a situation in which each individ- 
ual, upon arrival at a certain location, fixes the time of 
next jump not only for itself, but also for all of its de- 
scendants, and have shown that this leads to clustering 
phenomena analogous to those in a random environment. 



What we have found is that correlated behaviors analo- 
gous to those induced by the spatial structure of the en- 
vironment could be generated through inheritance, with- 
out the need to invoke any additional spatial interaction 
among individuals. This similarity between correlated 
behaviors is lost instead if an "antisocial" condition is 
considered, in which each individual migrates indepen- 
dently, following an internal clock. 

There is therefore a prescription problem in the way 
newly generated individuals are made to jump, that is 
absent in the normal diffusion case. It is important to 
notice that this implies differences in the transport prop- 
erties of the population, that show up already at the 
level of one-individual statistics. We mention that sim- 
ilar prescription problems were observed in [13], which 
dealt with a reaction - anomalous diffusion system, with 
reaction of coagulation type: A + A — >• A. 

We shall consider a simple, spatially homogeneous pop- 
ulation model, with Markovian births and deaths, and a 
migration dynamics of CTRW type [9, 14]. Following [3], 
we shall refer to the individuals in the simplified pop- 
ulation model as "bugs". A single offspring will be as- 
sumed in each birth event, with equal birth and death 
rates T for stationarity. Both superdiffusive and subd- 
iffusive regimes will be considered, although, as it will 
soon appear, the degree of sociality of the bugs affects 
the long-time population dynamics much more than the 
regime of anomalous diffusion. 



II. THE MODEL 

We briefly recall here the main properties of the CTRW 
model (more details in Appendix A). 

The model is a generalization of the random walk, in 
which the walker migrates through a sequence of inde- 
pendent jumps, separated by randomly chosen waiting 
times. The joint PDF (probability density function) for 
the jump length x and waiting time t can be taken in the 
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following form, in one dimension (D = 1): 

At^ / x 2 \ 

^ l) ~ t^KxT) exp I - 2{AxW 2 ) ' (1) 

where < p < 1 and < t fixes a microscale for the 
problem. In the absence of birth and death processes, 
both subdiffusive and superdiffusive behaviors may be 
generated, setting appropriately the jump length Ax(t) 
[9]. Focusing on the case of a power-law dependence of 
the typical jump length on the waiting time: 

Ax(t) ~ Ax{Q){t/At) v ' 2 , t » At, (2) 

we have basically the following possibilities, provided we 
set t ^> At and again neglect the effect of deaths [9, 11]: 

v<p^ (\x(t) - x{0)\ 2 ) ~ V, 

p<i/=> (\x(t) - x(0)\ 2 ) ~n v t v , (3) 

where we have introduced the generalized diffusivity 
n a = [Ax(0)] 2 / At a . We see that subdiffusion and su- 
perdiffusion are obtained respectively for v < 1 and 
v > 1. However, as described in detail in Appendix 
A, the presence of demography prevents the long-time 
asymptotics described in Eq. (3) to be ever reached. 
Thus, the long-time evolution of a population that was 
concentrated initially at a single spot, is described, in the 
presence of demography, by a normal diffusion process. 

Of course, with or without demography, normal diffu- 
sion can be equally generated by choosing a waiting time 
PDF that is not fat-tailed. 

We shall study the CTRW model with births and 
deaths, given an initial condition at time t = 0, in 
which the bugs are uniformly distributed in the domain: 
Pi(x, 0) = no, where p\ is the bug density. As the param- 
eters in Eq. (1) are chosen independent of the coordinate, 
the density will remain uniform also for t > 0; in par- 
ticular, imposing equality of the death and birth rates: 
pi{x,t) = pi{x,0) = uq. Demography however, through 
fluctuations, induces a secular component in the higher 
order correlations, that is a manifestation of clustering. 

To study the phenomenon, following [15], we will con- 
sider the two-bug density 

P2(xi, 2 , ti, 2 ) = P2(x, t) = nl + p 2c (x, t), (4) 

and we shall derive evolution equations for the connected 
component /02c(x, t) in the two opposite regimes of an- 
tisocial and inheritance governed, strongly social bugs. 
Notice that, contrary to [15], in which it was sufficient 
to consider the evolution of the one-time PDF /^(x; t, t), 
the non-Markovian nature of the CTRW will force us to 
study the full two-time statistical problem. 

III. THE ANTISOCIAL CASE 

In this case, the jump PDF is reinitialized at each birth 
event. The dynamics is that of reactions in the pres- 
ence of anomalous diffusion, that has been extensively 
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FIG. 1: Contribution to the 1- and two-bug densities by dif- 
ferent classes of family trees. Vertical lines ending with the 
small circle indicate a factor <I>. Letters A, B, and D indicate 
arrival, birth and death events. A line with a kink ending in 
A indicates the sequence of a waiting time and a jump, and is 
associated with a factor tp, as indicated in Eq. (12). A birth 
event is associated with a factor T. 



studied in the mean field regime (see e.g. [16] and refer- 
ences therein). Extension to two-bug statistics requires 
derivation of an evolution equation for the two-bug den- 
sity p 2 (x,t) = P2(^l,2,*l,2): 

p 2 (x,t) = / dn / dr 2 J? 2 (x,t -t)$(ti)$(t 2 ), (5) 
Jo Jo 

where 77 2 (x, t — t) is essentially the rate of formation 
(birth plus arrival) of bug pairs in (x, t — t) and the 
^(t"i,2) are the probabilities that the bugs arriving or 
born in (x, t — r) survive without jumping until t\$. 
Since e _rt is the survival probability in a time t, <&(i) = 
<&(i)e rt is the no-jump probability; from Eq. (1): 

= J At J Ax ip{x,r) - {At/ty. (6) 

Equation (5) is the generalization of similar equations 
connecting the one-bug density p\ and arrival/birth rate 
771 derived in [9] and [16] (see Appendix A). 

The bug density p 2 (x,t) can be represented as a sum 
of contributions by family trees, of the kind illustrated 
in Figs. 1 and 2. The small circles on top of the graphs 
identify the argument of the density function (the one- 
bug density p\ in the first line of Fig. 1 and p 2 in the 
second) and the vertical lines ending in the small circles 
identify the factors <l. The part of the graphs that is 
obtained chopping out these top lines accounts therefore 
for the rates 771,2- 

Clustering is accounted for by the connected part of 
P2 ■ The family trees responsible for this contribution are 
those summing to case (b) in Fig. 2. In this picture, 
graphs (c) and (d) in Fig. 2 appear to play an important 
role. It turns out that they act as a discrete source in Eq. 
(5), which we will see is associated with the fact that any 
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FIG. 2: Family trees contributing to the two-bug density 
P2(x, t); (a): disconnected part; (6): connected part p2c\ (c) 
and (c?) [in (d) we have ti > tH > ti]: family trees acting as 
discrete source for p2 C . Notice that chopping down branches 
in (6), we end up with family trees in either form (c) or (d). 



family tree in case (b) of Fig. 2 can be obtained adding 
or extending branches in a family tree in (c) or (d). 

The discrete nature of graph (c) can be understood, 
based on the direct parenthood relation implied between 
bugs 1 and 2, which is an intrinsically discrete property 
[in case (d), this direct relation exists only when t\ = 
ti] . These considerations prompt us to separate from the 
start discrete and continuous contributions r\i = rj2 lsc + 
r]2° nt , with 772 lsc accounting for graphs (c) and (d) in Fig. 
2, and to rewrite the corresponding contribution in Eq. 
(5) as a source term: 



s(x,t) 



/ dT i / ' 

Jo Jo 



dr 2 ^*' c (x,t-T)$(ri)$(75). (7) 



Let us evaluate explicitly this component. Let us start 
by evaluating the increment AN V that the family trees 
represented in case c of Fig. 2, produce on Ny — 
J v dxipi(x, t), with V some interval containing x\. The 
width Ly of the interval V will play the role of a coarse 
graining length for p 2 . 

To obtain AN V , we must calculate the product of the 
number of bugs generated in (xi,ti), surviving without 
jumps to ti, and the number of bugs originated after t\ 
from the parent bug (PB) (plus the PB itself, if still alive) 
that are counted in V. The number of bugs generated 
by each PB in (n, t\ + dri), that survive at x\ without 
jumps till t\, is rdri$(ti — n). Let us indicate with 
Gv(ti,T\,r) the number of bugs originating from this 
same PB after n (again, including the PB itself, if still 
alive), that are counted in V at time ti (r indicates the 
time of birth or arrival at x\ of the PB). Taking into 
account the rate 1)\{x\,t) of birth and arrival of bugs 
in (xi,t) and the probability $(n — r) of their being 
available for parenthood in {x\,T\), we obtain: 



an v = r 



f 1 f T1 

/ dn / dr T?i(a:i,T)&(ii -n) 
Jo Jo 

x *(Ti-r)G v (t 2 ,Ti,T). (8) 



Now, for iy large enough, all the descendants of the 
PB will remain in V, and, provided the birth and 



death rates are equal: Gy (£2 > t\ , r) — > 1 . This causes 
the two integrals in Eq. (8) to decouple. We have 
/ Tl dr t)i{xi,t)$(ti - r) = pi(a;i,ri) = n , the uni- 
form one-bug density [this is just the evolution equation 
for the one-bug density; see Eq. (Al)]. From here we get 
the expression: 

AN v (t) = n r f 1 dnMh-n). 
Jo 

We proceed in identical way to calculate the contribution 
AA^y from graphs (d) in Fig. 2. An analogous calculation 
leads to the result 

AN V = r / 1 dr^ f 2 dr m (xi, t)*(*i - r 2 ) 
Jo Jo 

x $(r 2 -r)Gy(i 2 ,ii,^). 

= n r / dr 2 $(t! - r 2 ). 
Jo 

From here we obtain for ANy = AN V + AN V , using Eq. 
(6): 



AiVy(t) = 2n r / dTi$(ti -n) ~n (rAi)'\ (9) 
Jo 

We see that most bugs contributing in to A7Vy(t) 

are born in an interval of width At before t\. Thus, for 
ti—t\ 3> At, the bugs in V will have dispersed (dying 
and reproducing) for a time ~ ti — t\, and for Eq. (9) 
to hold, we need Ly ^> X(t 2 — t\) where X{t) is the 
typical separation at time t of bugs originating from a 
common ancestor at time zero. We thus see that, if we 
are interested in time scales t 2 — t\ ^> At, and we are 
able to coarse grain pi at a scale Ly ;§> X(ti — ti), we 
can write for the source in Eq. (7): 



s(x,t) ~ no(rAt)"5(a;i -x 2 ). 



(10) 



Let us pass to analysis of r^ ™*, that can be decom- 
posed into contributions from births (B) and arrivals (A): 



cont 
72 



(11) 

with r\i B (x\,tx\xi,ti) = ri2 A (xi,ti;xi,ti). The com- 
ponents r/2 A ', . . . obey equations that are generalizations 
of those for one- bug [9, 16] (see Appendix A). It is con- 
venient to isolate from the start the connected contri- 
butions to rj AA ,r]2 B and r/.f' 4 , meaning that the contri- 
bution from the initial condition p 2 (x, 0) is disregarded 
(contribution by graph a in Fig. 2). This leads to the 
result: 

^(x.t) = / dn / dr 2 / d 2 y 77 2c (x-y,t-r) 
Jo Jo J 

X ^(j/l,n)^(j/2,T 2 ), 

^c 8 ^'*) = y dr y dy ?7i c (xi - - r|x 2 t 2 ) 

x ip(y,T), 

Vi c B (x,t) - r 2 p 2c (x,t), (12) 
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where ip(x,t) = t/j(x,t)e~ rt is the PDF that the bug 
makes a jump of length x at time t after birth or after 
a previous jump, and that therefore it did not die in the 
meanwhile. One more equation is needed for r]i c (xt\x' t') , 
that is the rate of birth and arrival of bugs at (x, t) given 
presence of a bug at (x', t'). This is the equation for the 
conditional density p-L C {%iti\x 2 t 2 ) = Uq p 2c (x, t): 

pic(x 1 t 1 \x2t 2 ) = [ dr r) lc (x 1 T\x 2 t 2 )$(ti - r). (13) 
Jo 

The integral equations (12) have a graphical interpre- 
tations as operations of adding (through birth) or ex- 
tending (through jumps) branches in the family trees of 
Figs. 1-2. For instance, the contribution ?72c( x — y,t — 
T)ip(yiTi)ip(y2 1 T2)d 2 Td 2 y to f7 2 c( x > t) can be seen as the 
operation of extending a branch from (x\ — yi,t\ — ti) to 
(xi,t\), and another one from (x 2 — y 2 ,t 2 — t 2 ) to (x 2 ,t 2 ) 
in the family trees contributing to 772c (x — y, t — t). This 
substantiates the role of graphs (c) and (d) in Fig. 2 as 
source for p 2 through Eqs. (5) and (12) [notice that 7/2 lsc 
still contributes in the RHS of the first and second of Eq. 
(12)]- 

A. Dependence on the kind of diffusion 

The system of equations (5,11-13) can be solved by 
Fourier-Laplace transform: /(x, t) — > /k z . Defining 
P2ckz = 27r<5(fci + k 2 )C klZ and s kz = 2ir5(ki + k 2 )a klZ , 
and indicating z lt2 = T + z\ t2 , the result is, after little 
algebra (see Appendix B): 

{1 - ipkz^kz^ - rbPkz 1 ®z 2 + i>kih®2i] 

+r 2 $ 5l $ 52 }C fcz = a kz . (14) 

We are going to consider separately the parameter ranges 
in Eq. (3) that would lead, in the absence of demography, 
to sub and super diffusion (the case of normal diffusion is 
going to be discussed in Appendix C). We thus consider 
H < 1, v = for subdiffusion, and p < 1, v > 1 for 
superdiffusion. Fourier-Laplace transforming Eq. (1) in 
the limit k, z\ t2 — > 0, we obtain the result, in the two 
cases: 

%l> kz = 1 - {MzY - (zAt)»- a (Axk) 2 /2, (15) 

where a = fi, v, respectively, in the subdiffusive and su- 
pcrdiffusive case, and we have reabsorbed in the two pa- 
rameters At ~ At, Ax ~ Ax(0) any numerical coefficient 
arising in the Fourier-Laplace transform. In the same 
way, we obtain from Eq. (6) the expression for <& z valid 
in the two cases: 

$ z = At***"- 1 . (16) 

If we are interested in evaluating p 2c (x 7 t) at separations 
t 2 — t\ ^> At and x 2 — x\ ^> X(t 2 — ti), we can use Eq. 
(10) to evaluate the source: Ofc z ~ n (rA£) M / (z\z 2 ); in 



the long-time regime Tt\ t2 3> 1, we can Taylor expand 
in z\ )2 , and Eq. (14) becomes in both superdiffusive and 
subdiffusive cases: 

C kz ~ p + Kfc 2 )z 1 z 2 ]- 1 n r, (17) 

where z — z\ + z 2 and R = KaT 1 ^", (again, a — ji, v in 
the subdiffusive and superdiffusive case). 

The reason why both subdiffusive and superdiffusive 
regimes lead to the same Eq. (17) is that in the limit 
rzi j2 — \ 0, we can expand $ Zl 2 ~ $r + &yZ\, 2 , and sim- 
ilarly for iphzi 2 - I n other words, the condition rzi.2 <S 1 
destroys any anomalous scaling with respect to z\^ 2 in 
2 and il>kzi,2, and therefore also in C kz . This has 
the counterpart in the original variables (x, t) in the ob- 
servation that for rti i2 » 1, the dominant factor in the 
individual bug dynamics is its own death, which prevents 
the long waiting times responsible for anomalous scaling 
in Eq. (3) to be ever reached (see Appendix A). 

Inverse Fourier-Laplace transform of Eq. (17) leads to 
the final result for C(x,t) = /?2c(x, t)\ X2 - Xl=x : 

C(x,t) ~ n Q X-\T^I 2 f{x 2 /{Kf)), (18) 

where i = (h + t 2 )/2, X 2 = RT' 1 - (^(r- 1 )) 2 is the 
square of the typical bug displacement in a lifetime, and 
f(x) is a function equal to one at x = and decaying 
rapidly to zero at \x\ > 1. The scaling in / shows that 
for \t 2 — ii| < ( it is possible to set Ly 3> X(\t 2 — t\\) 
and still have Ly much smaller than the spatial scale of 
variation of C at time i. By continuity, this allows us to 
evaluate C(x, t) at x = 0,ii = t 2 , where the derivation 
leading to Eq. (10) and then to (18) would not apply. 

We notice the physical meaning of the parameter R, 
that is the only place in which memory of the initial 
scaling exponents of the process is preserved. This is 
the effective diffusivity for the bugs in the presence of 
demography. Equation (18) is basically telling us that for 
Tt 3>, the bugs behave as if they were simple Brownian 
walkers, which execute at intervals T -1 steps of order 
X r . 

The scaling in Eq. (18) is the same as obtained in 
[15], in the case of individuals undergoing normal diffu- 
sion (Brownian bugs). As illustrated in Figs. 3 and 4, 
this result is confirmed by numerical simulations. The 
similarity with the results in [15] extends to D > 1. Mul- 
tiplying the right hand side (RHS) of Eq. (17) by k°- x 
and then inverse Fourier-Laplace transforming, we con- 
firm in fact the result in [15], of logarithmic divergence 
of C(x,t) for t — > 00 for D = 2 and no divergence for 
D = 3. 



IV. THE SOCIAL CASE 

The analysis carried on so far illustrates the impor- 
tance of initialization upon birth of the jumping PDF. 
One expects that anomalous behaviors may occur if a 
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FIG. 3: Result of numerical simulations in D = 1 for fj, = 
1/2, v = (without demography, this would correspond to a 
subdiffusive regime). Main graph: (a) plot of C — r~ (1+fl)/2 C 
vs. t = t/At in the antisocial case and fit by t 1 ^ 2 ; (s) plot of 
C = T~ C vs. t = t/At; in the social case; (d) plot of (7 = 
r _1 C vs. t — t/At; in the case of transport by a random field 
with correlation length l c = Ax/IG. The fit of (s) and (d) is 
provided by t. In all three cases (a), (s) and (rf), dotted, heavy 
and thin lines correspond respectively to TAt = 0.01, 0.05, 0.1. 
Insert: plot of the correlation length A (units Ax) of C in 
the antisocial case and fit by A oc t 1 ^ 2 . In all cases: total 
population iVo = 2 x 10 5 ; domain size L — 2 x W 4 Ax. In case 
(d) we have put e = 10 [see Eq. (31)]. 




10 2 10 3 10 4 10 5 10 6 
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FIG. 4: Same as Fig. 3 in the case fi = 1/2, v = 3/2 (with- 
out demography, this would correspond to a superdiffusive 
regime) . 



less strong form of reinitialization for tp is adopted. How- 
ever, it turns out that also imposing an excessive social 
constraint could result in anomalous scaling break-up. 

We consider a situation in which the time of first jump 
for a bug is determined by the ancestor that initiated the 
colony to which the bug belongs. All the bugs in that 
colony (and, if still alive, the ancestor itself) will jump 
therefore at the same time, independently. It is interest- 
ing to notice that this social regime is indistinguishable, 
at the one-bug level, from what would be obtained in the 
absence of demography. The reason for this is that the 
number of individuals in a colony - if the birth and death 
rates are equal - will remain equal to one on the aver- 



age. The dispersal of the bugs, when their jump time ar- 
rives, will produce therefore the same result, essentially, 
as what would be obtained from considering different re- 
alizations of the initial bug in the absence of births and 
deaths. (The same result would obviously be obtained 
considering the still tighter constrain that all the bugs 
jump together as a single entity). 

This situation leads to a strongly simplified popula- 
tion dynamics at the two-bug level, as compared to the 
antisocial case. We can write in fact, setting t 2 > t\\ 



p 2 (x,t) = / dn / dr 2 ^ A (x,t-r) 

x (F( tl -T))(F(t 2 -r)) 

x $(ti-Ti)$(t 2 -75) + a(x,t), (19) 

s(x,t) = 5(x!-x 2 ) / dr t]i(xi,t) 
Jo 

x (F(h - r)F{t 2 - T))$(t 2 - r), (20) 

where r/j 4 is the rate of arrival of individual bugs, and F(t) 
is the fluctuating number of descendants at time t of a 
bug alive at t = 0. The variable F describes a Galton- 
Watson process [17], which, for equal birth and death 
rates, is characterized by a constant average (F(t)) = 1, 
and a linearly growing correlation: 



(F{t!)F(t 2 )) = l + 2rmin(ti,t 2 ). 



(21) 



Notice the absence in Eqs. (19-20) of hats on the 3>'s, 
and of contributions AB,BA and BB in r\ 2 , as all de- 
mography is contained in the factors F. Absorption of 
demography in the factors F and the fact that (F) = 1 is 
what makes the dynamics identical at the one-bug level, 
to the one of a simple CTRW in the absence of birth and 
death. At the two-bug level of Eq. (19), however, we still 
find a source s(x, t), that acounts for the fluctuations in 
the number of individuals in a given colony. As illus- 
trated in Eq. (20), this quantity must be proportional to 
the arrival rate of individual bugs at x\ = x 2 , and to the 
probability that the colony does not disperse before the 
final time t 2 . 

Let us analyze in detail the structure of this source 
term and start by evaluating its long-time behavior. No- 
tice that, thanks to the fact that s(x, t) depends only on 
space difference, we can write s(x, t) — ak(t)6(xi — x 2 ). 

Let us focus first on the equal time case and proceed 
by Laplace transforming <Jk(t,t). We obtain: 

C z [a k ] = 2Tt 1 ^ , z ~ n rz- 2 => a k (t,t) - n Q Tt, (22) 

where use has been made of Eqs. (16,20) and (21), $' z = 
d§ z /dz, and of the result, generalizing Eq. (19) to one- 
bug statistics: n = J * dri]f(x, r)$(t — r) =>■ & z ilfz = 
n /z [compare with Eq. (Al)]. On the other hand we 
have, from Eqs. (20) and (21), for t 2 > t\: <r fe (t) — 
r$(t 2 ) L 1 drr]i (r)(ti — t). Laplace transforming in t\, 
using again Eqs. (16,20-21) and r\y z = no/(z$ z ), we get: 



C Zl [a k ] - n (] Tz 1 2 



CT fc (t) ~ n rt(ii/t 2 )". (23) 
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We could interpolate Eqs. (22-23) to obtain, for t 2 >t\: 

s(x,t) ~ n Tt 1 (t 1 /t 2 TS(x 1 - x 2 ). (24) 

We notice at this point the important difference with 
respect to the antisocial regime. While in that case fluc- 
tuation build-up was produced as a balance between a 
constant 0((At)^) forcing and diffusion, in the present 
case, this build-up is realized directly in s(x, t), which is 
here a macroscopic quantity [18]. In the present case, the 
only possible balance in Eq. (19) is p 2 — s at x x = x 2 , 
which tells us that p 2 (x, t) must grow linearly in time for 
x\ = x 2 . 

To be convinced of this result, let us write down ex- 
plicitly the evolution equation for p 2 . We separate in the 
arrival rate r\ 2 A the connected component rj^, which 
obeys an equation analogous to the first of Eq. (12): 

r?4 A (x,t) = / dn / dr 2 / d 2 y r? 2c (x - y, t - r) 
Jo Jo J 

x 4'{yi,T 1 )tp{y 2 ,T 2 ), (25) 

with rj 2c = r]^ + , q 2 hsc , and r(% ,sc providing the discrete 
contribution to the social dynamics. This last contribu- 
tion is associated with pairs of bugs that belong to a same 
colony and can be written in the form 

r# sc (x,t) = 2TS(t 1 ~t 2 ) J\t Jdyr)?(y,T) 

x ip(xi - y\ti - r)ip(x 2 - y\t 2 - t), 

x {h -r)V(ii -t), (26) 

where J?i(y,r) gives the rate of arrival of the common 
ancestors at (y,r), and 2T(ti — t), from Eq. (21), is 
the average number of pairs of descendants per ancestor, 
available to jump at time t\ = t 2 . [Of course, ip{t) is the 
waiting time PDF and ip(x\t) = i>( x >t)/' l l } (t) is the jump 
length conditional PDF]. 

Substituting Eqs. (25) and (26) into Eq. (19), ex- 
changing the order of the convolutions and setting again 
(F) = 1, wc obtain the expression: 

P2c(x,t) = / dn / dr 2 / d 2 y V(yi,n)-0(y2,T 2 ) 
Jo Jo J 

x K(x-y,t -t) + s(x,t), 
where K = p 2c + s — s, and 

s(x, t) = f 1 dn f 2 dr 2 r/f sc (x, t - t) 
Jo Jo 

x $(n)$(n). (27) 

From here we reach the final result: 

i o 2c (x, t) = s(x, t) + p 2c (x, t), 
where p 2c obeys the equation 



and 



(28) 



P2c(x,t) 



ti rt2 

dn / dr 2 / d z yip(y 1 ,Ti)ip(y 2 ,T 2 ) 
o Jo 



x p 2c (x-y,t-r) + s(x,t), 



(29) 



/•II /"I2 [■ 

s( x ,t) = / dr i / dr 2 / d 2 y ip(yi,T 1 )ip(y 2 ,T 2 ) 
Jo Jo J 

x B (x-y,t-T) (30) 

is a source term that, as can be checked from Eqs. (26,27) 
and (30), is continuous with respect to x\ — x 2 . 

Comparing with Eqs. (24) and (29), we thus see that 
the two terms on the RHS of Eq. (28) account respec- 
tively for the singular and regular components of p 2c . 
The singular component s(x, t), in particular, coincides 
with what would be observed in a population that does 
not migrate (in fact, its dynamics is that of a Galton- 
Watson process). The only place in the singular compo- 
nent, where migration seems to have some effect, is the 
structure of time correlations, as accounted for by the 
power-law decay in Eq. (24) [compare with Eq. (21)]. 
As illustrated in Figs. 3 and 4, the clustering dynamics 
at long times is dominated by the singular component 
s(x, t) in both the superdiffusive and the subdiffusive 
case. 

We stress that these results are conditioned to the wait- 
ing time PDF ip{t) being fat-tailed. If this condition were 
not satisfied, it is possible to see that p 2c would not have 
a singular component and the dynamics would become 
that of Brownian bugs [18]. 



V. POPULATION DYNAMICS IN A RANDOM 
ENVIRONMENT 

The kind of social behavior discussed in the previous 
section is very close to what one could expect if the bug 
migration were determined by a spatial arrangement of 
random traps. The opening time of the trap determines 
the time of dispersal of the bugs that are located in it. 
When a trap opens, all the bugs jump independently and 
remain caught either in the same trap (whose opening 
time is now updated) or in some other trap elsewhere. 

The important difference, with respect to the case 
of inheritance produced social behaviors, is that now, 
colonies that were initiated by bugs that happened to be 
caught in the same trap, will share an identical dispersal 
time. In the case considered in the previous section, the 
dispersal time of colonies sharing the same location in 
space was in general different. 

To verify the possibility of discrepancies between the 
two mechanisms for social behavior, a series of numerical 
simulation have been carried on in D = 1 and D = 2, 
with random traps whose opening times are indepen- 
dently distributed with PDF ip(t) oc t" 1+ ^. The jump 
length PDF for the bugs ip(x\t) = tp(x } t)/ip(t) is the one 
provided by Eq. (1). 

In all cases considered, periodic boundary conditions 
were utilized, and the total number of bugs was kept 
fixed, varying from the outside the birth and death rates 
in the population. (The meaning of a constant popula- 
tion constraint is to subtract any global finite population 



7 




ii 




y. ^ .". * - V*. 
■ ■■■ ■ • 

if: 


*; . '-»•■- 


"■/■*>"'■"■ 


v-: 














■^s.*'- : >*- 














£ ■ .>.•>'-:■ 





■ . •■ , i, 

■ --t-r , ■ •:. } ** 

. ■- '• V. «' 



FIG. 5: Sequence of snapshots of the population spatial dis- 
tribution from a numerical simulation of bugs in a D — 2 
random environment. Going clockwise from top left: t/At = 
4 x 10 3 ,4 x 10 4 ,4 x 10 5 ,4 x 10 6 . Parameters of the simu- 
lation: jj, = 1/2, v — 3/2 (would be superdiffusive case); 
TAt = 0.01; K = 10 _3 (Aa;(0)) 2 /At. Domain size in units 
Ax(0): 1400 x 1400. Correlation length for random field: 
l c = 14Aa;(0). Total population: N = 2 x 10 5 ; occupation 
number: e = 10. 



effect from the local fluctuations associated with cluster- 
ing [19]). In addition to this, in D = 2, the possibility 
of a small normal diffusivity k acting beside the jumps, 
has been considered, to mimick the effect of individual 
small scale motion of the bugs. (In the case of plankton, 
this additional component may be interpreted as due to 
individual swimming, while the long jumps are associ- 
ated with eddy break-up and and transport of plankton 
patches away, along strain filaments). A picture of what 
would occur in this case is presented in Fig. 5. 

It is clear that the discrepancy between different so- 
cial mechanisms will be negligible if the population is 
so dilute that on the average, at most one colony (and 
therefore one bug) is present at a given time in a given 
trap. Given a typical size l c of the trap, an important 
parameter is therefore the product 

e = n l? (31) 

that gives the typical occupation number of a trap. 

What we have found, in both D = 1 and D = 2, is 
that clustering is remarkably insensitive to whether the 
population is concentrated or diluted, or to the presence 
of small scale motions that cause bug leakages from the 
traps. As illustrated in Fig. 6, a change of more than two 
orders of magnitude in e does not produce any apprecia- 
ble change in the fluctuation amplitude C(0; t, t). In the 
same way, as can be seen in Figs. 3 and 4, the difference 
in the profile for C(0;t, t) in the two cases in which so- 
ciality is produced by inheritance and spatially extended 
random traps, are negligible. The same result holds both 




t 



FIG. 6: Plot of C = r _1 C vs. t = t/At and fit by t, in the 
case of transport by a random field, (a): D — 2, the same 
case as Fig. 5; the fluctuation amplitude C(0; t, t) was coarse 
grained at scale l c . The case Ko = leads to a profile that 
is indistinguishable from the one in figure, (b) D = 1; again 
a would be superdiffusive case, with fi — 1/2, v — 3/2 and 
TAt = 0.01; the three curves correspond to to e = 1 (heavy 
line), e = 20 (thin line), e — 200 (dotted line); for the rest, 
the same parameter as case (d) in Figs. 3 and 4 were utilized. 
A slope oc i is shown for comparison. 



in the subdiffusive and the superdiffusive case. 

An explanation for the insensitivity of clustering to the 
value of the parameter e, and to the mechanism for so- 
ciality, is the fact the clustering behavior accounted for 
by Eqs. (24) and (28) is already extremal. The maximal 
fluctuation growth, corresponding to individuals that do 
not migrate, as described by the Galton- Watson regime 
of Eq. (21), is already achieved in the inheritance gov- 
erned dynamics of the previous section. Passing to a 
random trap governed dynamics, further diminishes the 
degree of dispersal of the bugs and maintains the clus- 
tering on the maximal level C(0; t, t) oc t. Increasing the 
value of the parameter e goes in the same direction and 
does not modify the result. 

Geometrically, as accounted for by Eqs. (24) and (28), 
this corresponds to a situation in which, unless some 
small scale migration is assumed, the typical cluster size 
is zero. Repetition of the D = 2 simulation of Fig. 5 
with kq = can be shown to lead, in fact, to a bug dis- 
tribution in the long time limit, that is a discrete set of 
isolated towers, as it would occur with bugs that do not 
migrate at all [3] . The same can be seen to occur also in 
the D = 1 simulations, in both the inheritance and ran- 
dom trap governed cases. This contrasts with the result 
in the antisocial case (see insert in Figs. 3 and 4), which 
describes cluster whose size A is a growing function of 
time. 



VI. CONCLUSIONS 

The main result of the present analysis is that social 
behaviors can play a strong role in the transport prop- 
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erties of a population, whose individuals migrate follow- 
ing a CTRW strategy. An interesting aspect, with po- 
tentially useful application, is that enforcing a simple 
inheritance-social constraint on the bugs, is sufficient to 
replicate the dynamics of a population in a random envi- 
ronment (say, an assembly of random traps). This means 
that it is not necessary to impose spatial interactions 
among the CTRW-bugs to mimick the common response 
of nearby individuals to the environment. Nevertheless, 
a social constraint is necessary to guarantee that the dy- 
namics of CTRW-bugs and that of individuals in a field of 
random traps be identical. Antisocial bugs, for instance, 
do not satisfy this property and their dynamics coincides 
instead with that of Brownian bugs. In particular, as in 
the Brownian case of [15], antisocial bugs are character- 
ized by reduced clustering in D > 1, with fluctuations 
that grow only logarithmically (D — 2) or saturate to a 
constant level (D = 3). 

The clustering behavior in the social case (and there- 
fore also in the case of a random environment) coincides 
with what would be observed in the case of bugs that 
do not migrate. The growth of the fluctuation amplitude 
is described by a Galton- Watson process independently 
of the dimensionality of the problem, and the size of the 
clusters is determined by eventual small scale motions 
of the individuals. (This is opposite to the behavior of 
antisocial and Brownian bugs, where the clusters are pro- 
duced by a competition between migration and demog- 
raphy). In consequence of this, no clustering reduction 
with increasing dimensionality occurs in the social case. 

Thus, in realistic situations of a population dispersed 
in a bi- or three-dimensional random environement, a 
transport dynamics governed by random traps with long 
permanence times, will lead to much stronger clustering 
behaviors than it would occur if the transport dynamics 
were of the normal diffusion type. This result is inde- 
pendent of the kind of transport, either superdiffusive or 
subdiffusive, that would result if no birth and death ef- 
fects were present. The only important constraint is that 
the waiting time PDF in the CTRW (or in the random 
traps) be fat-tailed. On the contrary, the strong role of 
sociality is lost in the case of Brownian bugs. Also, trans- 
port by an assembly of random traps with exponentially 
distributed opening times, turns out to be described by 
a Brownian bug dynamics, without the need of imposing 
social constraints of any sort. 

In all cases, demography appears to be very good at de- 
stroying anomalous scaling behaviors. Antisocial behav- 
iors bring the dynamics back to that of Brownian bugs; 
strongly social behaviors bring it back to that of non- 
migrating individuals. In the antisocial case, anomalous 
scaling could possibly be recovered for t <C r _1 . From 
Eq. (18), we see however that fluctuations become rele- 
vant, i.e. C(x, t, t) ~ rip, only when Ti ~ (noX-p) 2 , which 
is small unless the density is so low that the typical bug 
separation is at least of the order of their displacement 
in a lifetime. 

The question remains open on whether anomalous scal- 



ing could be obtained from regimes half-way between the 
strongly social and antisocial cases that have been con- 
sidered so far. Anomalous scaling could arise as well from 
social behaviors that do not originate from inheritance, 
but are the result of more general spatial interactions be- 
tween the bugs, as discussed e.g. in [20]. Similarly, the 
question remains open on the effect of demography on 
bugs in which anomalously diffusive behaviors are pro- 
duced by means different from a CTRW, say, a fractional 
Brownian motion [21] or a sequence of Levy flights [20]. 
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Appendix A: Anomalous diffusion and demography 
in the antisocial case 

To check for the presence of anomalous diffusion be- 
haviors, we study the scaling of the bug density pi(x,t), 
given the initial condition pi(x, 0) = S(x). This choice of 
initial conditions makes pi(x,t) coincide with the PDF 
of finding the original bug or one of its descendants at 
position x at time t. 

The integral equation satisfied by pi has a form anal- 
ogous to Eq. (5) [9, 16]: 

Pl (x,t) = f dT Vl ( X ,T)^(t-T), (Al) 
Jo 

where r]\ is the rate of arrival and birth of individual bugs 
at (x, t), and $(£) = c~ rt <I>(i) is the probability that the 
bug does not jump or die in a time t [see Eq. (6)]. Let 
us decompose as in Eq. (11), rji into contributions from 
arrival and birth: r/i = + . The arrival contribution 
rji(x, t) is obained summing over all the bugs that where 
born or arrived at a previous time r at any position y 
and from there jumped at x at time t. Following [9], it is 
possible to include in (x, t) an initial condition p{x, 0), 
so that we can write, for t > 0: 

Vi(x,t) = J dr f dyrn(y,T)i>(x-y,t-T) 

+ P(x,0)6(t), (A2) 

while r](x,t) = for t < 0. The factor i>(x,t) = 
e~ rt ip(x 7 t) is the PDF that the bug jumps at at a time 
t after a previous arrival or birth, that the length of the 
jump is x, and that the bug is still alive at the time of 
the jump. The birth rate r]f(x,t) satisfies instead: 

» ? f(a;,t)=rp 1 (a;,t). (A3) 

Substituting Eqs. (A3) and (A2) into Eq. (Al), we ob- 
tain, after inverting the order of the convolutions, the 
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evolution equation: 

pi(x,t) = J dr J dy p l {y,r)ip{x - y,t - r) 

+ r Pl (x,T)$(t-T)\ +p(x,0)&(t). (A4) 

Fourier transforming in space and Laplace trasforming in 
time, Eq. (A4) becomes, imposing the initial condition 
pi(x,0) = S(x): 



Plkz = 



i - - 



where, as in Eq. (14), z = T + z. From here we obtain, 
in the the regime zAt -c 1, using Eqs. (15) and (16): 



Plkz 



1 1 

— 77 Kiy Z 

z 2 



l-a,-2 fe 2 



(A5) 



where n a ~ Ai 2 At _Q , with a = /x, f, respectively, in 
the subdiffusive and superdiffusive case. The variance 
X 2 (t) of the PDF pi is obtained as the inverse Laplace 
transform of the quantity 



d 2 pikz 



8k 2 



k=0 



c.r 1 -", r < z < Af- 1 



z < F 

where r 1 We obtain the final result 

'n a t a , Af<t<r- 1 



(A6) 



X\t) 



Kt, 



t > r- 1 



(A7) 



and we see that demography kills anomalous scaling for 

rt > i. 



Appendix B: Equations for the two-bug statistics in 
the antisocial case 

Let us fill the steps to go from Eqs. (5,11-13) to Eq. 
(14). Substituting Eqs. (12) into the connected part of 
Eq. (5), we find: 

pt\ />t2 p 

p 2c (x,t) = / dn / dr 2 d 2 y p 2c {y,r) 
Jo Jo J 

x ^(x! - j/i, ti - T 1 )tp(x 2 - yiM - T 2 ) 



+ 



Jo Jo 



dr 2 



[r/ 2 r(x,r)+^(x,r 



%T(x,r) $(ti-ri)$(i 2 -r 2 ) 
s(x,t), 



(Bl) 



where use has been made of Eq. (7) to eliminate rj c 
from the expression. On the other hand, we can write, 



using the second of Eq. (12) and Eq. (13): 

1 dn / 2 dr 2 ^ c B (x,T)$(ti - Ti)$(t 2 - r 2 ) 
o Jo 

Z"* 1 /" Tl /" 
= n r / dr 2 / dn / dr[ / dt/i 6(i 2 - r 2 ) 
Jo Jo Jo J 

x?7i c (yiTi|x 2 T 2 )-0(xi -yi,n - r()$(n -n) 



= r / dr 2 / dn / dyi $(t 2 - n) 

Jo Jo J 
xp2(yiTi;x 2 T2)^(xi - 2/1, n - n). 



(B2) 



Substituting Eq. (B2) with the third of Eq. (12) into 
Eq. (Bl), we get the result, after inverting the order of 
the convolutions: 

p 2c (x,t) = J dn J dr 2 | J d 2 y p 2c (y,r) 

x ■4>(x 1 - ?/i, n - Ti)i>(x 2 - y 2 M - r 2 ) 
+ r / dt/i p 2 c(yi,x 2 ;r)iil)(x 1 - yi,h -n) 

x $(i 2 - r 2 ) + r / dy 2 p 2 (a;i,2/ 2 ; r) 



x -0(x 2 - y 2 ,t 2 - r 2 )$(n - n) 
+ r 2 i o 2 (x,r)<l(n-Ti)$(t2-T 2 )} 

+ *(x,t), (B3) 

which, Fourier transforming in space and Laplace trans- 
forming in time, gives Eq. (14). 



Appendix C: Application to the case of Brownian 
bugs 

Using Eqs. (14) and (10) in the Brownian bug case is 
like using a shotgun to take down a mosquito. Neverthe- 
less, it allows one to test the techniques derived in the 
present paper on the known results in [15]. 

In the case of Brownian bugs we have, in place of Eqs. 
(1) and (6): 



^ l) = ^aJa x CXP ( - A~t 2A^) (C1) 



t x 2 \ 



and 



$(t) = exp(-t/At). (C2) 
From here we can write, in place of Eqs. (15) and (16): 

At 



C M-(Axk) 2 /2) , 

= — TTaTz — ' 



1 + Atz 



■ (C3) 



Substituting Eq. (C2) into Eq. (9), we obtain, in place 
ofEq. (10): 



s(x, t) = 2n TAtS(xi — x 2 ). 



(C4) 
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Substituting Eqs. (C3) and (C4) into Eq. (14) and Tay- 
lor expanding around k, z = 0, we find in place of Eq. 
(17): 

C fez = 2[{z + Kk 2 )z 1 z 2 }- 1 n T, (C5) 

where z — z\ + z 2 , and k = Ax 2 / At is the diffusion of 
the Brownian walkers. Contrary to the case of Eq. (17), 
the transport part of the dynamics, accounted for by the 
diffusivity k, is independent of demography We see that 
Eq. (C5) is the Laplace transform of the forced heat 
diffusion equation 

d t C(x, t) = nd 2 x C{x, t) + 2n a T8(x), (C6) 

where C(x, t) = C(x; t, t). 

It is possible to see that Eq. (C6) corresponds to the 
continuous limit of the result in [15]. Rather than tak- 
ing the continuous limit of that result, we take the limit 
directly on the original problem on the lattice. 

Let the discretization in time be given by At and 
the one in space by £ > Ax. Thus C(x,t) — 
l~ 2 (N(0, t)N(x, t)), with N(x, t) the instantaneous num- 
ber of bugs in the cell x. We can easily fix the jump 
probability, to obtain a diffusive dynamics with diffusiv- 
ity k = Ax 2 /At. Analysis of the source term in Eq. (14) 



requires a little more work. Let us indicate by AC(x,t) 
the contribution to C(x,t + At) — C(x,t) from demog- 
raphy. Again, we easily see that for x ^ 0, due to in- 
dependence of the birth and death processes, we have 
£ 2 AC(x,t) = (N(0,t)AN(x,t)) + (AN(0,t)N(x,t)) = 
0. If x = 0, however, we have that the quantity 
(AA 2 (0, t)\ N(0, t)) is non-zero. Introducing the indica- 
tor function Si = 0, 1, 2 depending on whether bug i dies, 
remains inactive, or gives birth, we have in fact a contri- 
bution to (AN 2 (0,t)\N(0,t)): 

N(0,t) 

(S i S j ) = 2N(0,t)TAt, 

ij 

where we have used the fact that (SiSj) = for j / j, 
and Si = 0,1,2 with probabilities TAt, 1 - 2FAt and 
TAt. Exploiting the fact that (N(0,t)) = not, we find in 
the end 

AC(0,i) = 2n Q TAte-\ 

which leads in the continuous limit to the expression for 
the source provided in Eq. (C6). 
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